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We present the first next-to-leading-order QCD analysis of neutrino charm production, 
using a sample of 6090 z/^- and 77^-induced opposite-sign dimuon events observed in the 
CCFR detector at the Fermilab Tevatron. We find that the nucleon strange quark content 
is suppressed with respect to the non-strange sea quarks by a factor k, = 0.477 
where the error includes statistical, systematic and QCD scale uncertainties. In contrast 
to previous leading order analyses, we find that the strange sea a;-dependence is similar 
to that of the non-strange sea, and that the measured charm quark mass, rric = 1.70 ± 
0.19 GeV/c^, is larger and consistent with that determined in other processes. Further 
analysis finds that the difference in ^-distributions between xs{x) and xs{x) is small. A 
measurement of the Cabibbo-Kobayashi-Maskawa matrix element \Vcd\ — 0.232 iao2o is 
also presented. 
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1. Introduction 

Nucleoli structure at high momentum transfers is characterized by parton distribution 
functions, which describe the proton and neutron in terms of quarks and gluons using the 
factorization theorems of Quantum Chromodynamics (QCD) These nucleon parton 
distributions are essential inputs when using perturbative calculations to predict high en- 
ergy processes involving nucleons, such as those at the Tevatron collider and the planned 
LHC. The ability to consistently predict such processes using one set of universal par- 
ton distributions is an important test of QCD as the theory of the strong interactions. 
Neutrino-nucleon deep-inelastic scattering is particularly suited for measuring the parton 
densities due to the neutrino's ability to resolve the flavor of the nucleon constituents. 
Furthermore, neutrino scattering is an effective way to study the dynamics of heavy quark 
production, due to the light to heavy quark transition at the charged current vertex. In 
particular, neutrino charm production can be used to isolate the nucleon strange quark 
distributions, xs{x) and xs{x). 

We present the first next-to-leading-order (NLO) QCD analysis of neutrino and an- 
tineutrino production of charmed quarks including the first direct determination of the 
strange quark distribution defined at NLO. The order of the strange quark distribution 
from this analysis matches that of recent global nucleon structure analyses, like those 
of CTEQ 0] and MRS [Q. Since the data presented here provide the most sensitive 
constraints on the strange quark distribution, these results should become important in- 
gredients in future global parton distribution fits. 

In addition, we present measurements of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix element \Vcd\ and the mass of the charm quark rric- The charm quark mass from 
this analysis is directly comparable with rric measurements derived from NLO analyses 
of other processes. Such comparisons are a test of the perturbative QCD phenomenology 
for heavy quarks. Neutrino production of charm off nucleon d quarks provides a direct 
determination of the product |V^rfp-Bc, where Be is the weighted average of the semi- 
leptonic branching ratios of the charmed hadrons produced. \Vcd\ is isolated by including 
an estimate of Be using data from other experiments, including a re-analysis of Fermilab 
neutrino-emulsion data. 

The signature for the production of charmed quarks in neutrino- and antineutrino- 
nucleon scattering is the presence of two oppositely-signed muons. In the case of neutrino 
scattering, the underlying process is a neutrino interacting with an s or c? quark, producing 
a charm quark that fragments into a charmed hadron. The charmed hadron's semileptonic 
decay produces a second muon of opposite sign from the first. 

i^/. + N — >/i-+c + X 
The analogous process with an incident antineutrino proceeds through an interaction with 
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an s or (i antiquark, again leading to oppositely-signed muons in the final state. 

17^ + N — ^/i++c + X 

2. The CCFR detector and event selection 

The dimuon data were accumulated during two runs, E744 and E770, with the Chicago- 
Columbia-Fermilab-Rochester (CCFR) detector at the Fermilab Tevatron Quad- Triplet 
neutrino beam. This wide-band beam was composed of and with energies up to 600 
GeV and a flatter energy spectrum than characteristic of horn-focused neutrino beams. 
In the CCFR detector 0, neutrino interactions occur in the 690 ton unmagnetized 
steel-scintillator target-calorimeter, which is instrumented with drift chambers for muon 
tracking. The calorimeter's hadronic energy resolution is a/E = 0.85/ \/E. The target 
is followed by a solid-iron toroidal magnetic spectrometer, which identifies muons and 
measures their momenta with a resolution Ap = O.llp. 

The detector measures and p^j, the momenta of the two muons, 6*^1 and 6*^2, 
the angles of the muons with respect to the neutrino beam axis, and E^^^, the en- 
ergy of the hadronic system. Muon 1 is labeled as the leading muon, defined as the 
one emerging from the leptonic vertex, using a procedure described below. The visi- 
ble energy in a dimuon event, E^i^ = E^-^ + E^^ + E^^^, misses the energy of the decay 
neutrino. Therefore, we make the distinction between the visible quantities — quantities 
derived directly from measurement — and the physical quantities, which are inferred on 
average by correcting distributions of the visible quantities using the Monte Carlo sim- 
ulation described below. Variables commonly used to describe deep inelastic scattering 
are: Ql^^ = 4£'vis-E'^iSin^(^^j/2), the negative square of the four-momentum transfer, 
a^vis = Qvis/[2-^(-^had + -^^2)]' ^he Bjorken scaling variable, where M is the nucleon mass, 
?/vis = (-Ehad + -E^2)/-^vis, the inelasticity, and = M'^ + Ql-^^{l/x^is - 1), the invariant 
mass squared of the hadronic system. 

In the experiment, charged-current single muon events are required to have i?vis > 30 
GeV, E'hid > 10 GeV, Ql^^ > 1 GeV^/c^ and p^^ > 9 GeV/c. Dimuon events are selected 
by making the further requirement that the second muon has > 5 GeV/ c and that both 
muons have < 0.250 rad. The second muon's momentum is measured in the magnetic 
spectrometer whenever possible, otherwise it is determined from the muon's range in the 
target. In order to reduce non-prompt sources of second muons, events in which muon 
2 does not reach the toroid must also satisfy -Ehad < 130 GeV. The final dimuon sample 
contains 6090 events and is characterized by (-Evis) = 192 GeV, (W^is) = 168 GeV^/c^, 
(Qlis) = 25.5 GeV^/c^, and (xvis) = 0.15. Results from a leading-order analysis of this 
data sample were reported previously 0. 
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3. Differential cross section 



The differential cross section for dimuon production is expressed generally as 

D{Z) 5e(c^/i+X), 



d^a{u^N ^ fi-fi+X) (Pa{v^N cX) 



dy dz d^ dy 

where the function D{z) describes the hadronization of charmed quarks and is the 
weighted average of the semi-leptonic branching ratios of the charmed hadrons produced 
in neutrino interactions. 

The heavy charm quark introduces an energy threshold in the charm production rate. 
This is a kinematic effect for which ^, the momentum fraction of the struck quark, is 
related to the Bjorken scaling variable x, through the expression |^ 




(2) 



where rric is the charm quark mass and iris refers to the initial state quark mass, either 
the strange quark or the down quark, and A = A(— Q^, 'm^, m^) is the triangle function, 

defined by A(a, 6, c) = ^Ja^ + + — 2{ab + bc + ca). The full expression for ^ can be 
simplified by neglecting the small effect of the initial state quark mass to yield 

i^4^^4M^-'^\- (3) 



Q 

Relating ^ and x through the charm quark mass is referred to as slow-rescaling [§] . 

At leading order (LO) charm is produced by scattering directly off of strange and down 
quarks in the nucleon. The LO differential cross section for an isoscalar target, neglecting 
target mass effects, is given by: 



(Pa{u^N^cX)\ _ G^ME, ^ ^^^^^^ ^ , cmc „2m it/ |2 

cd 



+ 2«.K.,^)|K.P}(l-^). (4) 



where ^u{^,ij?), id{^^iJ?) and C,s{^,fi'^) represent the momentum distributions of the 
u, d and s quarks within the proton (the corresponding z7^ process has the quarks re- 
placed by their antiquark partners) and \Vcd\ and \Vcs\ are the CKM matrix elements. 
The dependence of the parton distributions on the scale fi^ is specified by QCD In 
the leading order analysis of Reference 0, Callan-Gross violation is included by replac- 
ing the term [1 - ml/{2ME^^)] in Equation (D with {[1 + Rl{^, /i"^)] [1 + {2M^/Q)Y^ 
[1 — y — Mxy/ {2E)] +xy/C,}, and using external measurements of the structure function 
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-Rl(C)/^^) HOI- NLO formalism, violation of the Callan-Gross relation emerges as 

a consequence of QCD. 

The LO expression illustrates the sensitivity of the process to the strange quark sea. 
Charm (anticharm) production from scattering off d (d) quarks is Cabibbo suppressed. In 
the case of charm produced by neutrinos, approximately 50% is due to scattering from s 
quarks, even though the d quark content of the proton is approximately ten times larger. 
In the case of antineutrino scattering, where d quarks from the sea contribute, roughly 
90% is due to scattering off s quarks. 

Because neutrino charm production has a large sea quark component at leading-order, 
the next-to- leading-order gluon-initiated contributions are significant |Tl|] . The size of the 
gluon distribution, which is an order of magnitude larger than the sea quark distribution, 
compensates for the extra power of involved in the gluon-initiated diagram. The NLO 
quark-initiated diagrams, shown in Figure Qd, in which a gluon is radiated, also enter the 
perturbative expansion at 0{as), but the contributions of these diagrams to the cross 
section are not enhanced by large underlying parton distributions. Calculations including 



the next-to-leading-order formalism have recently become available |T^, [T^, |1J] and lead 
to the analysis in this Letter. 





Figure 1: Mechanisms that contribute to 
neutrino production of charm up to 
a) The dominant diagrams: the leading- 
order quark-initiated diagram, and the t 
channel and u channel gluon-initiated dia- 
grams, respectively, b) The radiative-gluon 
and self-energy diagrams. 
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4. Monte Carlo simulation 



Information about the strange sea, the charm quark mass, and the branching ratio is 
extracted by comparing the x^is and E^is distributions of the data to theoretical expec- 
tations contained in a Monte Carlo simulation. The Monte Carlo program models the 
dependence of these physics parameters as well as the effects of detector acceptance, reso- 
lution smearing, and missing energy associated with the charmed particle decay. Dimuon 
Monte Carlo event generation proceeds by using a simulated sample of charged-current 
single muon events and demanding that the hadronic system contains a second muon from 
charm decay. The single muon Monte Carlo sample is normalized to the charged-current 
data sample, ensuring that the Monte Carlo energy spectrum exactly models that of the 
data. 

The species of charmed particles produced in neutrino interactions as a function of 



neutrino energy was measured by Fermilab E531 |T5|. With an E^, > 30 GeV cut, the 
production is dominated by charged and neutral D mesons. Fragmentation to D^s in 
the Monte Carlo simulation is parameterized by the fragmentation function of Collins 
and Spiller H, D{z) = [(1 - z)/z + e(2 - z)/{l - z)] (1 + zfll - (l/z) - e/(l - z)]-^, 
where z = pn/pf)^^ is the fraction of its maximum momentum that the D meson carries 
and e is a free parameter. The parameter e in the fragmentation model is fit by using the 
distribution of ^^is = E^^/{E^^ + Ehad)- 

The dimuon events are divided into those from incident or by a separation 
procedure that assumes that the leading muon has larger transverse momentum with 
respect to the direction of the hadron shower than the muon from the charmed hadron 
decay. This procedure separates the sample into 5030 z/^-induced events and 1060 V^- 
induced events. The largest uncertainty in this procedure is due to the knowledge of 
the charm meson p± distribution. For this study we parameterize the p± distribution by 
dn/dp^ cx e~^^i and use the Fermilab E531 emulsion data JlSl with W"^ > 30 GeV^/c^ 



to determine f3 = 1.21 ± 0.34. Using this method, the separation procedure is found to 
misidentify 5.8 ± 0.4% of the u events and 7.3 ± 0.4% of the u events. 

The charm-initiated dimuon signal is contaminated by non-prompt pion and kaon 
decay. The high-density calorimeter minimizes this contamination due to the short inter- 
action length of the detector. A combination of hadronic test beam muoproduction data 
and Monte Carlo simulations predicts a small vr/K decay background of 797 ± 118 z/^ and 
118 ±25 events 0. 

To calculate the probability of producing charm, we employ the NLO QCD charm 
production differential cross section calculation of Aivazis, Collins, Olness and Tung 
ll2| , including the Born and gluon-fusion diagrams, shown in Figure ||a. These are 



the leading contributions to charm production. The calculation is performed in the MS 
scheme. The factorization scale in the calculation is chosen to be /i = 2pf^^, where 
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^max _ ml, M^) / a/ AW^ is the maximum available transverse momentum of the 

initial state quark coming from the gluon splitting, or equivalently of the final state charm 
quark, for the given kinematic variables x and Q^. The renormalization scale is chosen to 
equal the factorization scale. We discuss the uncertainty due to the choice of these scales 
below. Electromegnetic radiative corrections to the cross section are calculated using the 



method of Bardin et al. 18 



The finite momentum cut on the second muon limits the acceptance of events at- 
tributable to gluon-initiated production, 

W^g cs. 

Gluon-initiated production of charm proceeds through both the t and u channels as shown 
in Figure |l|a. While these diagrams are quantum mechanically equivalent, they dominate 
different regions of phase space. In the t channel, the gluon splits into an ss pair and 
the c quark emerges from the W-boson vertex. In the u channel, the legs of the c and 
s quarks are crossed — the gluon splits into a cc pair and the s quark emerges from the 
W-boson vertex. 

In the W-boson-gluon center of mass frame, the c quark is produced at an angle 6* 
relative to the W-boson direction. The production angle is related to the momentum of 
the c quark in the lab, and hence with p^j. When 6* is small — t channel dominance — 
the c quark carries most of the W-boson momentum. As 6* approaches vr — u channel 
dominance — the c quark emerges with little momentum in the lab. Consequently, events 
with large 6* are less likely to produce a second muon with > 5 GeV/c. 

An acceptance correction due to this effect is determined by folding the calculated 
squared production amplitude for producing charm at angle 9* with the experimental 
acceptance. To determine the experimental acceptance, the ability of events with finite 6* 
to pass the 5 GeV/ c cut is compared to that when the events are generated with 9* = 0. 
We find that this relative acceptance drops to near zero at 6* = ir, and is about 74% at 
9* = tt/2. Integrating over the calculated squared production amplitude, which peaks 
in the forward and backward directions, the overall t-u channel acceptance correction 
for gluon-initiated production is 60 ± 10%. The effect of this acceptance correction is 
small but not insignificant; for example, it shifts the value of rric determined from the fit 
described below by +0.07 GeV/c^. 

Measurements of the F2 and xF^, structure functions by CCFR [|T^, are used to 
determine the singlet and the non-singlet quark distributions, xqsi{x,iJ?) = /i^) -|- 
xq{x, /i^) and xqNs{x, /i^) = xq{x, iJ?) — xq{x, /i^), respectively, and the gluon distribution, 
xg^x,^"^) ||2T|. These distributions are obtained from next-to-leading-order QCD fits to 
the structure function data using the QCD evolution programs of Duke and Owens 

Pi- 
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To resolve the strange component of the quark sea, the singlet and non-singlet quark 
distributions are separated by flavor. Insofar as isospin is a good symmetry, our experi- 
ment is insensitive to the exact form of the up and down valence and sea quark distribu- 
tions, because the neutrino target is nearly isoscalar. An isoscalar correction accounts for 
the 5.67% neutron excess in the target. 



The proton valence quark content, xqv{x,ii^) = xqj\fs{x, fi'^), is parameterized by 

a;gy(x,yU^) = xuyix, fx^) + xdv{x, jj^), 
xdv{x, jji^) = Aci{l — x)xuv{x, fx^), 



(5) 



where the shape difference for xdv{x) better fits charged- lepton scattering measurements 
of F2/F2 [^. is fixed by demanding that the ratio of the number of c? to m valence 
quarks in the proton is 1/2. 

The non-strange quark and antiquark components of the assumed to be sym- 

metric, so that xu^x^jj?) = xus{x,fi'^), = xds{x,iJi?). The isoscalar correction 

is applied assuming xu^x^jj?) = xd^x^jj?). The strange component of the quark sea is 
allowed to have a different magnitude and shape from the non-strange component. The 
strange quark content is set by the parameter 



/q^[xs(x, /i^) + xs{x, ^'^)]dx 
/q^ [xu{x, /i^) -|- xd{x, /i^)] dx ' 



(6) 



where k = 1 would indicate a flavor SU(3) symmetric sea. The shape of the strange 
quark distribution relates to that of the non-strange sea by a shape parameter a, where 
a = would indicate that the strange sea has the same x dependence as the non-strange 
component of the quark sea. Shape parameters are defined for each of the two fits below. 
In the first fit the strange quark and antiquark distributions are assumed to be the same; 
in the second fit a;s(a:, /i^) and /x^) are fit separately. 

4.1 xs{x,fi'^) = xs{x,n'^) fit 

This fit assumes that xs{x, /i^) and xs{x, /i^) are the same. The sea quark distributions 
are parameterized by: 



xq[x^ /i 



xs{x, n^) 



xu{x, 11^) + xd{x, jj? 



X] 



2 

xu{x, /i 



+ xs{x, fi 



'^^ + Xd{x, yU^) 



(7) 



where A^ is defined in terms of n and a. 



A minimization is performed to find the strange sea parameters k and a, the values 
of Be and rric, and the fragmentation parameter e, by fitting to the a;vis, E^is and Zvis 
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distributions of the dimuon data, projections of which are shown in Figure |[ Taking 
iKdl = 0.221 ± 0.003 and = 0.9743 ± 0.0008 12D as input values and using the 
Colhns-Spiller fragmentation function, the extracted NLO parameters with their statis- 
tical and systematic errors are presented in the first line of Table |1|. In this primary fit, 
the uncertainties due to fragmentation are included in the statistical errors of all of the 
parameters. The value of = 52.2 for 65 degrees of freedom suggests excellent agreement 
between the data and the NLO theoretical model. 

Our previous LO results 0, which were found by fitting to the Xvis and i?vis dis- 
tributions of the same data sample and using the Peterson fragmentation function ||26|| , 
D{z) = N{z[l - (l/z) - ep/(l - z)]'^}-^ with ep = 0.20, are listed in the third line 
of Table 1. For comparison with these results. Table 1 includes the NLO parameters 
determined using the same fit procedure. The two fits using the Peterson fragmentation 
function include the uncertainty due to fragmentation in the systematic errors of all of 
the parameters. 





Fragmentation 




K 


a 


Be 


me (GeV/c^) 


NLO fit 


Collins-Spiller 


52.2/65 


0.477 


-0.02 


0.1091 


1.70 




e = 0.81 ±0.14 


+0.046 +0.023 
-0.044 -0.024 


+0.60 +0.28 
-0.54 -0.26 


+0.0082 +0.0063 
-0.0074 -0.0051 


, Q .7 +0.09 
^^■'^1 -0.08 


NLO fit 


Peterson 


41.2/46 


0.468 


-0.05 


0.1047 


1.69 




ep = 0.20 ±0.04 


+0.061 +0.024 
-0.046 -0.025 


+0.46 +0.28 
-0.47 -0.26 


±0 0076 +0-0065 
itU.UUm _o.0052 


±0 16 


LO fit 


Peterson 


42.5/46 


0.373 


2.50 


0.1050 


1.31 


Ref. [6] 


ep = 0.20 ±0.04 


1HS± 0.018 


+0.60 +0.36 
-0.55 -0.25 


±0.007 ± 0.005 


+0.20 +0.12 
-0.22 -0.11 



Table 1: Next-to-leading-order and leading-order fit results, assuming xs{x) = xs{x). Er- 
rors are statistical and systematic, except that the errors on the fragmentation parameters 
are statistical only. 



Estimates of the systematic uncertainties are obtained by varying model parameters 
within errors and are itemized in Table 0. 

Parton distributions are defined to a given order and scheme in QCD. Therefore, 
the magnitude of a given parton distribution differs between leading-order and next-to- 
leading-order. At NLO, the nucleon strange quark content is found to be k = 0.477 Ig^oso) 
indicating that the sea is not SU(3) symmetric — qualitatively the same result as from the 
LO analysis. The strange quark content may alternatively be given by 



J^[xs{x,iJ,'^) + xs{x, iJ,'^)]dx 
Jq[xu{x, fi'^) + xd{x, fi'^)]dx' 



so that by comparing to the total non-strange quark content, rj is less sensitive to changes 
in the determination of the sea quark content alone. At /i^ = 22.2 GeV^/c^ the ratio of 
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Figure 2: Xvis, -E-vis and z^is distributions for u- and 77-induced dimuon events. Data 
are given by the points and the sohd histogram is the result of fitting the dimuon event 
simulation. The dotted histogram is the background contribution to the former from pion 
and kaon decay. 
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Ul 


-De 


lll'C 


/I / XV UdL/JY^i U LliiU. 


+ 0.0010 
- 0.0022 


+ 0.150 
- 0.160 


- 0.0031 
+ 0.0027 


+ 0.006 
- 0.003 




- 0.0051 
+ 0.0078 


- 0.048 
+ 0.059 


+ 0.0028 
- 0.0002 


+ 0.031 
- 0.057 


Tfi 1 Q "1" "WTC^ ll rTrCl"t'l/~\Tl 

1 tJldLl Vc UdllUi dLlUli 


- 0.0059 
+ 0.0020 


+ 0.115 
- 0.077 


+ 0.0013 
- 0.0002 


+ 0.031 
- 0.007 


LiC LjCL/ ulwli Ciiit^ldil^y _l L /o 


- 0.0070 
+ 0.0021 


+ 0.044 
- 0.003 


+ 0.0022 
- 0.0012 


- 0.007 
+ 0.013 




+ 0.0014 
- 0.0032 


- 0.042 
+ 0.006 


+ 0.0009 
- 0.0014 


+ 0.036 
- 0.030 


F2 and xFs 


+ 0.0143 
- 0.0068 


+ 0.084 
- 0.114 


+ 0.0024 
- 0.0019 


+ 0.053 
- 0.017 


xqsiix), xqNsix) and 


+ 0.0110 
- 0.0110 


+ 0.180 
- 0.180 


+ 0.0009 
- 0.0009 


+ 0.038 
- 0.038 


iKsl and \Vcd\ 


+ 0.0095 
- 0.0099 


- 0.024 
+ 0.031 


- 0.0028 
+ 0.0032 


- 0.011 
+ 0.014 


t-u channel accep. corr. 


- 0.0135 
+ 0.0052 


- 0.009 
+ 0.062 


- 0.0002 
+ 0.0010 


- 0.020 
+ 0.008 



Table 2: Sources of systematic uncertainty in the determination of the fit parameters. 



antiquarks to quarks in the nucleon at NLO is found to be / dxxq{x, ji^) / J dxxq{x, fi"^) = 
Q/Q = 0.245 ± 0.005 and thereby the strange quark content with respect to the non- 
strange quarks is 

r/ = 0.099 ±0.008 ±0.004 ;°:°°^. (9) 

Since a nonzero value of a would indicate a shape difference between xq{x) and xs{x), 
the value a = —0.02 1 o;go indicates no shape difference at NLO. At leading order, we find 
the strange quarks softer than the overall quark sea by a factor (1 — x)° with a = 2.5±0.7. 
The difference in a between NLO and LO is attributable to the NLO xq{x) being softer 
than its LO counterpart, as shown in Figure ^. Figure ^ shows the NLO and LO xs{x), 
again indicating that the NLO distribution is larger in magnitude and softer than its LO 
counterpart. 

The strange quark distribution from the xs{x, /i^) = xs{x, jj?) fit is tabulated in Table 
^ for a few values of x and /x^ and is plotted at = 4 GeV^/c^ in Figure ^. The 
distribution can be parameterized by a function of the form a (1 — xY' x~'^. The values of 
the coefficients a, h and c are tabulated in Table ^. 

The charm quark mass parameter from the NLO fit is 1.70 ± 0.19 GeV/c^, which 
differs from the leading-order result, indicating the marked dependence of rric on the 
order to which the analysis is done. The NLO value of mc can be more consistently 
compared with measurements derived from other processes involving similar higher-order 
perturbative QCD calculations. A photon-gluon-fusion analysis of photoproduction data 
finds m, = 1.74tg:J| [p. 



The values of from the NLO and LO fits are consistent, providing a good check 
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(GeVVc2) 


X 


xs{x, /x^) 


1.0 


0.01 


0.126 ±0.012 ±0.006 


+ 0.008 
- 0.004 




0.05 


0.097 ±0.008 ±0.004 


+ 0.006 
- 0.003 




0.10 


0.068 ± 0.005 ± 0.003 


+ 0.004 
- 0.002 




0.20 


0.032 ± 0.003 ± 0.001 


+ 0.002 
- 0.001 


4.0 


0.01 


0.178 ±0.016 ±0.008 


+ 0.011 
- 0.005 




0.05 


0.111 ±0.009 ±0.005 


+ 0.007 
- 0.003 




0.10 


0.072 ±0.005 ±0.003 


+ 0.004 
- 0.002 




0.20 


0.030 ±0.003 ±0.001 


+ 0.002 
- 0.001 


20.0 


0.01 


0.229 ±0.020 ±0.010 


+ 0.014 
- 0.007 




0.05 


0.122 ±0.010 ±0.005 


+ 0.007 
- 0.004 




0.10 


0.073 ± 0.006 ± 0.003 


+ 0.004 
- 0.002 




0.20 


0.028 ±0.003 ±0.001 


+ 0.002 
- 0.001 


100.0 


0.01 


0.271 ±0.024 ±0.012 


+ 0.017 
- 0.008 




0.05 


0.128 ±0.010 ±0.005 


+ 0.008 
- 0.004 




0.10 


0.072 ±0.006 ±0.003 


+ 0.004 
- 0.002 




0.20 


0.026 ±0.003 ±0.001 


+ 0.002 
- 0.001 



Table 3: Values of xs(a;, /x^), defined at NLO using the MS renormalization scheme, from 
the xs{x,ii^) — xs{x,ii^) fit. The first error is statistical, the second is experimental 
systematic and the third is due to QCD scale uncertainty. 

of the fit procedure and the physics model. This value is constrained by neutrino charm 
production at x > 0.3 where the valence d quark contribution dominates, and thus is in- 
dependent of the strange and other sea quark distributions. The value of B^. = 0.109 1q;oo7 
agrees with an indirect determination, Bl — 0.099 ± 0.012, which is described in section 
5. 

As with all apphcations of perturbative QCD, a theoretical uncertainty is associated 
with the choice of factorization and renormalization scales. Some scale dependence is 
unavoidable for any calculation done to finite order in as- The /j,^ scale is interpreted as 
setting the boundary between the collinear and noncollinear regions of the p± integration 
over the final states. Therefore, a scale proportional to pY"^ is suggested by the authors of 
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Figure 3: The quark sea distribution Figure 4: The strange quark distribu- 
xq^x^jji^ = 4.0 GeV^/c^) determined at tion xs(x, /z^ = 4.0 GeV^/c^) determined at 
next-to-leading order and leading order. next-to-leading order (described in section 

4.1) and leading order. The band around 
the NLO curve indicates the ±lo" uncer- 
tainty in the distribution. 



/i^ GeVVc2 


a 


b 


c 


1.0 


0.135 


6.48 


0.000 


10.0 


0.117 


6.67 


0.096 


20.0 


0.107 


6.98 


0.164 


100.0 


0.100 


7.28 


0.210 



Table 4: The coefficients a, b and c from the parameterization of xs{x, /i^) using the form 
a(l — x)''x~'^, as described in the text. 



Ref. [IT], Figure ^ shows the scale dependence of the differential cross section, where 
the abscissa is in units of p™'^^. The scale dependence is weak for /i values above one unit 
of p™^, and there is a stronger scale dependence when /i is below this value. We choose 
/i = 2^™*^^ and find the scale uncertainty by varying /i between p™'^^ and 3p™^^. Fit results 
with various choices of the common factorization and renormalization scale are presented 
in Table |^. The values for these fits — all with 65 degrees of freedom — indicate that 
the data favor /x^ scales with smaller magnitudes. It should be noted that the values of 
the fit parameters are fairly insensitive to the choice of scale. 
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choice of scale, //^ 






a 


Be 


rUc (GeV/c2) 




50.4 


0.513 


0.18 


0.0987 


1.71 




52.2 


0.477 


-0.02 


0.1091 


1.70 




54.4 


0.460 


-0.10 


0.1142 


1.68 




51.7 


0.423 


-0.37 


0.1074 


1.80 




56.1 


0.410 


-0.46 


0.1159 


1.71 




59.4 


0.408 


-0.54 


0.1206 


1.73 




52.5 


0.421 


-0.03 


0.1066 


1.65 


1 iQ' + n^^) 


57.1 


0.109 


-0.16 


0.1151 


1.64 


+ (2me)2 


52.8 


0.428 


0.00 


0.1068 


1.62 


4 + (2me)2] 


57.3 


0.415 


-0.15 


0.1161 


1.63 



Table 5: Central values of the fit parameters for various choices of the QCD scale 
Each fit contains 65 degrees of freedom. 
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b 

t3 



X 



b 



0.5 



0.0 



0.5 



0.0 



y=0.5, E=200 GeV 



=0.5, E=200 GeV 



TJ 
x=0.05 



x=0.25 



x=0.15 



Figure 5: The //^ scale dependence of the 
differential cross section for neutrino and 
antineutrino production of charm, where /i^ 
identifies the factorization and renormaliza- 
tion scales. The scale /j, on the abscissa is 
in units of p^'"''- For E = 200 GeV and 
y = 0.5, the x = 0.05, 0.15, 0.25 lines cor- 
respond to p^^"^ = 6.6, 6.2, 5.8 GeV/c, re- 
spectively. 



4-2 xs{x) 7^ xs{x) fit 

Theoretical work has explored the possibility that the nucleon contains a sizable heavy 
quark component at moderate x — the possibility of so-called intrinsic heavy quark states 
within the nucleon [^]. Postulating intrinsic strange quark states leads to the prediction 
that the s quark momentum distribution will be harder than the s quark distribution ||30|| . 
We explore this possibility by performing a fit in which the momentum distributions of 
the s and s quarks are allowed to be different. For this study the sea quark distributions 
are parameterized by: 

n. ^ ( xuix, u^) + xdix, u^)\ xs(x, ii'^) + xs(x, a"^) 
xq{x,n) = 2 + , 
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xu{x, 11^) + xd{x, 11^) 



xs{x,ii'^) = As{l-xY 

xs{x,ii^) = A'^{l-xY' 
The s and s are constrained to have the same number 



xu{x, IJ?) + xd{x, 11^) 



[ s{x,n'^)dx= [ s(x,fj?)dx. 
Jo Jo 



(10) 



:iii 



As and A'^ are defined in terms of k, a and a'. 

In order to reduce the number of free parameters, this fit constrains the average 
charmed hadron branching ratio to the value obtained from other measurements, Bl — 
0.099 ± 0.012 (see Section 5.). We fit for four parameters: the strange quark parameters 
K, a, and Aa = a — a' and the charm quark mass rric. The result is: 



K 

a 
Aa 



0.536 ± 0.030 ± 0.036 ; ± 0.009, 
-0.78 ± 0.40 ± 0.56 ± 0.98 ± 0.50, 
-0.46 ± 0.42 ± 0.36 ± 0.65 ± 0.17, 

± 0.01 GeV/c^ 



1.66 ±0.16 ±0.071°:°^ 



(12) 



where the first error is statistical, the second is systematic, the third is due to the uncer- 
tainty in Sj, and the fourth is the error due to //^ scale uncertainty. 

The value of Aa = —0.46 ± 0.85 ± 0.17 indicates that the momentum distributions 
of s and s are consistent and the difference in the two distributions is hmited to —1.9 < 
Aa < 1.0 at the 90% confidence level. This is the first quantitative comparison of the 
components of the s and s quark sea. 

We also checked the assumption that the same average semileptonic branching ratio 
applies to the u- and 77-induced samples. A two parameter fit finds the branching ratio of 
i/-induced events = 0.1147 ± 0.0056, and AB^ = B^- B'^ = 0.011 ± 0.011, where B'^ 
is the branching ratio for I7-induced events and the errors are statistical only. The result 
indicates that there is no significant difference in the semileptonic decays of charmed 
particles and antiparticles at these energies. 

5. \Vcd\ measurement 

If the CKM matrix elements are not assumed, then the four parameter NLO fit in 
section 4.1 is performed by fitting a, rric and the following products: 



Wr 



K 



K + 2' 



cd\ Be 



\Vr.<,\ Br 



+ 0.38 
0.39 



= (5.34 
= (2.00 ±0.10 



0.27 
0.21 

+ 0.07 
- 0.05 



0.25^ 
0.5lJ 

0.06 
0.14 



X lO"'^ 
) X 10" 



(13) 
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These combinations can be used to extract [V^dP and fi:|V^sp when is determined 
from other data. Be is determined by combining the charmed particle semileptonic branch- 
ing ratios measured at e~^e~ coUiders ||2^ with the neutrino-production fractions measured 



by the Fermilab E531 neutrino-emulsion experiment [|1^]. Using an E'vis > 30 GeV cut, 
E531 determined the following production fractions: 52±6% 42±6% D+, 1±2% D+, 
and 5 ± 3% A+. In the E531 analysis, events that could not be unambiguously identified 
as or were all categorized as events. To remove this small bias, a re-analysis 
was performed that included updated values of the charmed hadron lifetimes This 
re-analysis finds the following production fractions with an i^vis > 30 GeV cut: 60 ± 6% 
26 ± 6% D'^, 7 ± 5% D'^, and 7 ± 4% A+. These production fractions are consistent 



with those measured by e+e experiments ||33| 



We find Bl = 0.099 ± 0.012 and extract the value of the CKM matrix element 

|Kd|= 0.232 (14) 

where the error indicates all sources of uncertainty, including the /i^ scale uncertainty. It 
compares very well with the PDG value, \Vcd\ = 0.221 ± 0.003, which is determined from 
measurements of the other matrix elements and the unitarity constraint on the CKM 
matrix assuming three generations. A measurement of \Vcs\ will be possible when an 
independent measurement of the strange sea content is available. 

6. Summary 

We have perfomed the first NLO QCD analysis of neutrino charm production and have 
measured the nucleon strange quark distribution and the electroweak parameters rUc and 
\Vcd\. We find: 

rUc = 1.70 ±0.19 ±0.02 

ft. — yj.^l I _ 0.050 + 0.036 

a = -0.02^°;^^^°;°^ (15) 

where k is the strange quark content with respect to the non-strange sea and a indicates 
a shape difference between the strange and non-strange sea, xs{x) oc (1 — x)"" [xu{x) + 
xd{x)]/2. The first error combines statistical and systematic errors in quadrature and the 
second is the uncertainty due to QCD fi"^ scale. This value of k indicates that the quark 
sea is not flavor SU(3) symmetric. The relative shape parameter a is consistent with 
zero, indicating that there is no shape difference between the strange and non-strange 
components of the sea. The value of rric obtained from the NLO analysis is consistent 
with that found in other processes. 

Using an externally determined production weighted charmed hadron branching ratio, 
Bl = 0.099 ± 0.012, we measure the CKM matrix element 

iKdl = 0.232 (16) 
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We have also studied the possibihty of a shape difference between the xs{x) and xs{x) 
distributions. We find that 



Aa = a-a' = -0.46 ±0.85 ±0.17, (17) 

where a and a' are shape parameters for xs{x) and xs{x), indicating no shape difference 
between the components of the strange quark sea. A shape difference is hmited to —1.9 < 
a — a' < 1.0 at 90% confidence level. 

The strange sea can also be inferred from a comparison of charged lepton and neutrino 
structure functions. To leading-order the lepton and neutrino structure functions are 
related by the "5/18ths rule," F!f / F^^ = 5/18{ 1 - 3/5[(s ± s)/{q + q)] }, where the 
strange sea enters as a correction. Comparison of structure function measurements from 
CCFR (z/Fe) P, H with those form SLAC (eD) H, NMC H and BCDMS (/iD) |35|, 



shows good agreement for x > 0.1 but a small discrepancy is seen between the neutrino 
and muon results for lower x. The source of this disagreement is under investigation. 



but may be due to the extra axial vector component present in neutrino scattering ||36|| . 
In contrast, the recent global fits by the CTEQ Collaboration (CTEQl distributions) 
have attributed the muon versus neutrino difference to an enhanced strange sea at low 
This possibility is ruled out by the measurements presented here. 

Acknowledgement 

We thank M.A.G. Aivazis, F.I. Olness and W.-K. Tung for encouragement, useful 
discussions and for providing computer code of their cross section calculation. We thank 
the management and staff of Fermilab, and acknowledge the help of many individuals 
at our home institutions. S.R.M. acknowledges the support of the Alfred P. Sloan and 
Cottrell Foundations. This research was supported by the National Science Foundation 
and the Department of Energy of the United States, who should be credited for their 
continued support of basic research. 



References 

[1] J.C. Collins, D.E. Soper and G. Sterman, in: Perturbative Quantum Chromodynam- 
ics, ed. A.H. Mueller (World Scientific, Singapore, 1989). 

[2] J. Botts et al. (CTEQ Collab.), Phys. Lett. B 304 (1993) 159. 

[3] A.D. Martin, R.G. Roberts and W.J. Stirling, Phys. Lett. B 306 (1993) 145. 

^The more recent CTEQ2 distributions include a preliminary version of the NLO strange sea results 
presented here, but do not account for the structure function discrepancy. 



19 



[4] W.K. Sakumoto et al. (CCFR Collab.), Nucl. Instrum. Methods Phys. Res. A 294 
(1990) 179. 

[5] B.J. King et al. (CCFR Collab.), Nucl. Instrum. Methods Phys. Res. A 302 (1991) 
254. 

[6] S.A. Rabinowitz et al (CCFR Collab.), Phys. Rev. Lett. 70 (1993) 134. 

[7] M.A.G. Aivazis, F.I. Olness and W.-K. Tung, MSU-HEP 93/15. 

[8] R.M. Barnett, Phys. Rev. Lett. 36 (1976) 1163. H. Georgi and H.D. Politzer, Phys. 
Rev. D 14 (1976) 1829. 

[9] V.N. Gribov and L.N. Lipatov, Sov. J. Nucl. Phys. 15 (1972) 438; G. AltareUi and 
G. Parisi, Nucl. Phys. B 126 (1977) 298. 

[10] S. Dasu et al, Phys. Rev. Lett. 61 (1988) 1061. 

[11] M.A.G. Aivazis, F.I. Olness and W.-K. Tung, Phys. Rev. Lett. 65 (1990) 2339. 

[12] M.A.G. Aivazis, J.C. Collins, F.I. Olness and W.-K. Tung, SMU-HEP/93-17. 

[13] G. Kramer and B. Lampe, Z. Phys. C 54 (1992) 139. 

[14] J.J. van der Bij and G.J. van Oldenborgh, Z. Phys. C 51 (1991) 477. 

[15] N. Ushida et al. (E531 Collab.), Phys. Lett. B 206 (1988) 375; B 206 (1988) 380. 

[16] P. Collins and T. Spiller, J. Phys. G 11 (1985) 1289. 

[17] P.H. Sandler et al. (CCFR Collab.), Z. Phys. C 57 (1993) 1. 

[18] D. Bardin and N. Shumeiko, Sov. J. Nucl. Phys. 29 (1979) 499. 

[19] P.Z. Quintas et al. (CCFR Collab.), Phys. Rev. Lett. 71 (1993) 1307. 

[20] W.C. Leung et al. (CCFR Collab.), Phys. Lett. B 317 (1993) 655. 

[21] The parton distributions are parameterized at a reference scale /Zg by xqNs{x, fi^) = 
Ansx'^'C^-x)'^^ +Bnsx'^'', xqsi{x,nl) = xqNs{x, nl) +Asea{l-x)'''''- 8indxg{x,fil) = 
Ag{l — x)^\ and evolved to all other values of //^ using the GLAP [^ evolution 
equations. Ajys is constrained by the fermion conservation sum rule / —xq^si^) = 3, 
and Ag is constrained by the momentum sum rule / dx [xqsi{x) + xg{x)] = 1. The 
final fit parameter is the strong coupling constant, given by Aj^, for a total of eight 
fit parameters. At /Xq = 1 GeV^/c^, the parameters are Aj^ = 0.210 ± 0.021 GeV/c, 
771 = 0.947 ± 0.027, T]2 = 3.81 ± 0.059, r/3 = 0.405 ± 0.047, Bns = 0.0210 ± 0.0039, 
Asea = 1.404 ± 0.047, r/sea = 6.49 ± 0.22, rjg = 3.18 ± 0.70, where the errors are 
statistical only. 



20 



[22] W.G. Seligman et al. (CCFR Collab.), in: Proc. XXVIIIth Rencontre de Moriond: 
QCD and High Energy Hadronic Interactions (Les Arcs, March 1993), ed. J. Tran 
Thanh Van (Editions Prontieres, Gif-sur-Yvette, 1993) p. 39. 

[23] D. Duke and J. Owens, Phys. Rev. D 30 (1984) 49. 

[24] P. Amaudruz et al, (NMC Collab.), Nucl. Phys. B 371 (1992) 3. 

[25] K. Hikasa et al. (Particle Data Group), Phys. Rev. D 45 (1992). 

[26] C. Peterson et al, Phys. Rev. D 27 (1983) 105. 

[27] J.C. Anjos et al (E691 Collab.), Phys. Rev. Lett. 65 (1990) 2503. 

[28] J.C. Collins, in: Proc. XXVth Rencontre de Moriond: High Energy Hadronic In- 
teractions (Les Arcs, March 1990), ed. J. Tran Thanh Van (Editions Prontieres, 
Gif-sur-Yvette, 1990) p. 123. 

[29] S.J. Brodsky, C. Peterson and N. Sakai, Phys. Rev. D 23 (1981) 2745. 

[30] M. Burkardt and B.J. Warr, Phys. Rev. D 45 (1992) 958. J.D. Bjorken, private 
communication. 

[31] T. Bolton, Nevis R#1501. 

[32] D. Bortoletto et al (CLEG Collab.), Phys. Rev. D 37 (1988) 1719. 

[33] A. Bodek et al, Phys. Rev. Lett. 50 (1983) 1431; 51 (1983) 534. 

[34] P. Amaudruz et al (NMC Collab.), Phys. Lett. B 295 (1992) 159. 

[35] A.C. Benvenuti et al (BCDMS Collab.), Phys. Lett. B 223 (1989) 485; B 237 (1990) 
592. 

[36] A. Donnachie and P.V. Landshoff, Z. Phys. C 61 (1994) 139. 



21 



